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^ ■ Abstract 

Ph ■ We consider the equations of relativistic magnetohydrodynamics (RMHD) in the case of 

0^ ■ special relativity. For the fluid rest frame a nonconservative reformulation of the RMHD 

equations gives a symmetric system for the vector of primitive (physical) variables. By 
applying the Lorentz transformation to this system we find a concrete form of symmetric 
matrices in the LAB-frame. The resulting symmetric system in terms of primitive variables is 
important for the study of various initial boundary value problems for the RMHD equations. 
We also find a so-called secondary symmetrization whose direct consequence is the extension 
of the sufficient stability condition obtained earlier for non-relativistic planar current-vortex 
sheets to the relativistic case. As in non-relativistic settings, this implies the local-in-time 
^ , existence of corresponding smooth nonplanar current- vortex sheets. 
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1 Introduction 



The equations of relativistic magnetohydrodynamics (RMHD) for an ideal fluid read (see, e.g., 

mm)- 

V„(pn") = 0, V„r"^ = 0, V„('u"6''-n'^6") = 0, (1) 

X: 

$H . where is the covariant derivative with respect to the Lorentzian metric g = diag (—1, 1, 1, 1) 

of the space-time with the components gajs (in this paper we restrict ourselves to the case of 
special relativity), p is the proper rest-mass density of the plasma, are components of the 
4- velocity, 

T"^ = {ph + 52)n"n^ + qg''^ - b^b^ , 

h = 1 + e + (p/p) is the relativistic specific enthalpy, p is the pressure, e = e{p, S) is the specific 
internal energy, S is the specific entropy, = b'^ba, b" are components of the magnetic field 
4- vector with respect to the plasma velocity, and q = p+ ^B^ is the total pressure. The 4- vectors 
satisfy the conditions u'^Ua = — 1 and u'^ba = 0, and the speed of the light is equal to unity. 

Let {x^,x) be inertial coordinates, with t = a time coordinate and x = {x^,x'^,x^) space 
coordinates. Then 

n° = -uo = r, u^ = Ui=Tvi, u := {ui,U2,U3) =Tv, 
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= 1 + 6° = -60 = (n, H), U = hi = ^ + (u, H)v^, 

b := {hiMM). = _ 5^ = M + Hf > 0, 

where F = (1 — lup)^-*^/^ is the Lorentz factor, v = (^1,^21^3) is the plasma 3-velocity, H = 
{Hi, H2, H3) is the magnetic field 3- vector, the Latin indices run from 1 to 3, and by ( , ) we 
denote the scalar product. 

The RMHD equations ([T]) can be rewritten as the system of conservation laws: 



dtipF) + div (pu) = 0, (2) 
dt {phru + \H\^v - {v,H)H) 

+ div{{ph + B'^)u(E)u-b®b) +Vq = 0, (3) 

dtiphr^ + -q)+ div {phFu + \H\'^v - {v, H)H) = 0, (4) 

dtH -V X (vxH) = 0, (5) 

where dt = d/dt, V = {81,82,83), di = 8/8x^, etc. Moreover, the first equation in the Maxwell 
equations containing in ([T|), 

divF = 0, (6) 
should be now considered as the divergence constraint on the initial data U{0,x) = Uq{x) for the 



unknown vector U = [p, u, H, S) which is one of the possible sets of so-called primitive (physical) 
variables [1]. 

Using ([6]) and the additional conservation law (entropy conservation) 

dt{pTS) + d\-v{pSu) = Q (7) 

which holds on smooth solutions of system ©-(IS]), and following Godunov's symmetrization 
procedure [6j , we can symmetrize the conservation laws ©-([Si) in terms of a vector of canonical 
variables Q = Q{U). This was done by Rugger i and Strumia [9] and also by Anile and Pennisi 
[21 [T]. However, a concrete form of symmetric matrices was not found in ^2]. Moreover, if we 
deal with an initial-boundary value problem, in particular, a free boundary problem it is very 
inconvenient to work in terms of the vector Q. 

In principle, from Q we can return to the vector of primitive variables U keeping the symme- 
try property (see jS])- But, for RMHD finding a concrete form of symmetric matrices associated 
to the vector Q and returning from Q io U \s connected with unimaginable (or even almost 
unrealizable in practice) calculations. To overcome this serious technical difficulty our main 
idea is to symmetrize the RMHD system for the fluid rest frame and then use the Lorentz 
transformation. Namely, taking into account ([6]), for the fluid rest frame {v = 0) we can rewrite 
©) ©, ([5]), in a nonconservative form which will be already a symmetric system. After 
that we should properly apply the Lorentz transformation to this system to get a corresponding 
symmetric system in the LAB-frame. 
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The main result of this paper can be formulated in the following theorem. 



Theorem 1.1 The RMHD equations ([2])- 
the symmetric system 

3 

Ao{U)dtU + J2^j(-^)^j^ = 
i=i 

for the vector U = {p, u, H, S), where 



on smooth solutions are equivalently rewritten as 

(8) 



An 



' ^ ^ 

pa'' 

V A 

0X0 

1 



A, 



ej 0\ 

ej Aj Mj'^ 

Mj VjM 

V v,J 



(9) 



°j = {Sij , 52j , Ssj) are the unit column vectors, c? = pp[p^S), 

A= iphT+^-^j I - (phr+^-^ ] V ^ V - - H (g) H + ^^^{v H + H v) , 



1 
f 



1 V ±± 

M = -{I + u®u), Mj = -b®ej--^h^v-^L 



El 

r2 



Aj = Vj\[ phT + ^ 1 / 



1^12 _ J^2\ 

phT + - — ^— \ v^v - -H ® H 

H ( 1 

+ Y2{v(S> H + H (g)v) -2{v,H){I -v(g>v) 



+ ^ {H ® Cj + ej H) - — {v (g) Cj + ej (g) v) 



and I is the unit matrix. The hyperbolicity condition Aq > holds provided that 

p> 0, Pp> 0, < < 1, 
where Cg = a'^/h = Pp/h is the relativistic speed of sound. 



(10) 



Clearly, we can use the symmetric form ([8]) for solutions with discontinuities (shock waves 
and characteristic discontinuities) in each smoothness domain as well as, for example, for the 
relativistic plasma- vacuum interface problem [13] , where the motion of the relativistic plasma is 
described by the RMHD equations. In fact, this interface problem can be mentioned as the first 
application of the symmetric form ([8]). The usage of ([8]) in [13] was really one of the key points 
in obtaining a priori estimates for the corresponding linearized initial boundary vale problem. 

In this paper we apply the symmetric form ([8]) to derive a sufficient stability condition for 
planar relativistic current-vortex sheets. As for the non-relativistic case [101 lllj. this stability 
condition being satisfied at each point of an initial nonplanar current-vortex sheet guarantees 
the local-in-time existence of current-vortex sheet solutions of the original nonlinear system. 
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However, for relativistic current-vortex sheets the symmetric form ([8]) is not enough for 
finding the sufficient stability condition and deriving a priori estimates for the hnearized problem. 
As in \10\ lllj. we need a so-called secondary symmetrization. Such a secondary symmetrization 
enables one to make the boundary conditions of the linearized problem dissipative. We call 
it secondary because the RMHD system being written in form ([SD is already symmetric, but 
there is a certain transformation of this system which keeps the symmetry property and has 
the advantage of dissipative boundary conditions. Actually, the hyperbolicity condition for 
the secondary symmetrization gives us the sufficient stability condition for planar relativistic 
current- vortex sheets. We write down it in the same form as in the non-relativistic case \10\ [TT] 
(see CTl. 

The plan of the rest of the paper is the following. In Sect. 2 we prove Theorem ll.il i.e., we get 
the desired symmetric form ([8]). In Sect. 3, we obtain another (secondary) symmetric form of the 
RMHD system with a more restrictive hyperbolicity condition. At last, in Sect. 4 we formulate 
the free boundary problem for relativistic current-vortex sheets, write down a corresponding 
linear stability problem for planar discontinuities, and deduce the above mentioned sufficient 
stability condition. 



2 Symmetrization of RMHD in terms of primitive variables 

Note that equations ([2]) and ([7]) imply 



dS 
'dt 



0, 



(11) 



with d/dt = dt + {v, V). In the absence of the magnetic field, H = 0, following [12], from ([2])-(l4 
(jlip we obtain the nonconservative form of the relativistic Euler equations 

r dp 
pa? dt 



2 j-i- + ^*^) + div u = 0, 



, , f du ( duW , „ dS 

Equations (I12p . ()13p form the symmetric system 



(12) 
(13) 



for W = {p, u, S), where the matrices Ba were written in |12j: 



Bn 



V phVB 




V 



B, 



pa'^ 



Cj phujB 



V 



V 







B 



(14) 



V QSi V. 



The natural idea of symmetrizing the RMHD equations in the general case of a nonzero 
magnetic field is the repetition of simple arguments above that could give the desired symmetric 
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system ([8]) for U = {p, u, H, S), and it is reasonable to expect that the matrices with cancelled 
fifth, sixth and seventh rows and columns coincide with in ()14p for H = 0. Unfortunately, 
this way is connected with extremely long technical calculations. To avoid them we propose the 
following procedure. First, it is enough to get the linearized counterpart 

3 

Ao{U)dt{5U) + Y,Aj{U)dj{5U) = (15) 
i=i 

of system ([8|), where 5U = {5p, 5u, 5H, 5S) is the vector of perturbations and U is now a constant 
vector, and p = p{p,S), e = e{p,S), Vi = Ui/T, F = (1 — w^)~^/^, 1/a^ = pp{p,S) and h = 
1 + e + (p/p) are constants. Below we write down system (|15p for the fluid rest frame {v = 0). 
After that we should only properly apply the Lorentz transformation and get system in the 
LAB-frame. 

The last equation in (jlSp will be the linearization of (jlip . 



dt{dS) + {v,V{6S))=0, (16) 
and it is natural to believe that the matrices An in (1151) will have the block structure 



Aq 0\ ^ A^ 0\ 

That is, our goal now is to find the symmetric matrices A^ of the system 

3 _ 

Ao{U)dt{5V) + Y,Aj{U)d,{5V) = (18) 
i=i 

for 5V = {5p,6u,5H). 

Let U' = {p,0,H',S') = U\^=o and 6V' = {6p',5u',6H') = 5VU=o. It is clear that the 
thermodynamical values coincide in the rest and LAB-frames: p' = p, p' = p, h' = h, 6p' = 5p, 
etc. We can easily write down the linearized system 

-^dt{6p') + dw{5u')=0, (19) 
pa^ 

{ph + H''^)dt{5u') - H'{H',dt{5u')) + V{5p') + H' X (V x {5H')) = 0, (20) 

dtiSH') - {H', V)((5n') + H'div {6u') = (21) 



associated with equations ([2]), ([3|) and ([5]) in the rest frame. Here we have taken into account (|16p 
for V = and the divergence constraint div {SH') = 0. Equations (|19p - (|2ip form the symmetric 
system 

_ 3 _ 

A',{U')dt{6V') + A'^{U')d,{6V') = 0, (22) 
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with 



1 



■ 

^' 
\00lJ 



fo e] 



A', 



\ 



e, JCj 
ICj 



(23) 



where 



A' = [ph + H'^) I-H' ]Cj = -H'jl + H' ® 

and the magnetic field is transformed under the Lorentz transformation as 

The last formula is obtained by applying the Lorentz transformation to the 4- vector {b^ , b) : 



(24) 



L 





^ r - 1 , ^, 



where 



L 



r -Tv 
r - 1 

-Tv I H t: — t; (g) f 



L- 



;) 

Tv 



. r- 1 

i H ;^ — f (Ki w 



Now we should properly apply the Lorentz transformation to system ()22|) . The perturbations 
are transformed as 

8V' = J5V (25) 
and it is clear that the matrix J has the form 

/l o\ 



J 



(26) 



Ji 

VO J3 J2/ 

where Ji is found by applying the Lorenz transformation to the perturbation ((JF, 6u) = {{v, 6u), 6u) 
of the 4- velocity (F,u), and J2 and J3 are found by the same procedure applied to the pertur- 
bation {6b^,6b) of the 4- vector {b^,b). Namely, we have: 

/ 



L 



L 



{v, 5u) 
6u 







1 - F 

, 5u + {v,5u) ^ „ V 



' 



(n, 6H) + {H, 6u) 
5b 



-{H,6u) 
, J3 6u + J26HJ 



5H' ) 



where 



Ji 



1 - F 



■vi^v, J2 = - ( / 



F - 1 



■ V (>$! V 



(27) 
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1 1 — r 

J3 = ®v + {v, H')I + —5- {v, H')v (g) v, 



and using (pll) . we recalculate J3: 



J3 = if u + (w, H)I - (u, (g) V. 



(28) 



In view of (I25p . 



and for calculating the elements c^^ of the matrices Ca through the elements a^^ of the matrices 
A'^ we apply the Lorentz transformation to the 4- vector {a\i^^\i^'^'ku^\i)' 



L 



-1 



VliJ 



That is, 



r - 1 



,2 ^j- 



1 - r 



( r?;jyi[, + A'j + — ^ vj-e ] J = vjAo + vj^^j'^gj + j^^;. J, 



(29) 
(30) 



where 





















3 


V 







, Co = r 


pa? 

V 


A' 




, 1C = H' (g)v - {v,H')I. 


k=l 




/c 


J 




\ 


K 


1 1 





Using we recalculate: 



IC = l;H (g)v - {v, H)I + ^— f (?;, F)r; (g) u, 
1 i f ^ 



r- 1 



1 1 - r 

ICj = --{Hjl - H ej) + -^-^ {v, H){vjl - v ej), 

A' = {ph + B'^) -l^H(^H + ^^{v,H){v(g)H + H(^v) +^^^^^{v,Hfv0v. 
Then, from (jl7p . ()23p . (|26p - ()30p after long calculations we find the symmetric matrices A^ in 

dH). 

As is known [IJ, natural physical restrictions guaranteeing the hyperbolicity of the RMHD 
system do not depend on the magnetic field and coincide with corresponding ones in relativistic 
gas dynamics. In our case, by direct calculations one can show that the hyperbolicity condition 
j4o > is equivalent to the condition Bq > for the relativistic Euler equations and holds 
if inequalities (jlOp are satisfied (of course, by default we also assume that \v\ < 1). The last 
inequality in p2p is the relativistic causality condition. 
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Remark 2.1 Strictly speaking, to prove the equivalence of system ([2]) -([5]) and ([8]) on smooth 
solutions we should also derive equations ([I])-® from system To this end we write down 
the subsystem for H contained in ([8]) and deduce from this subsystem the divergence constraint 
([6|) provided that it was satisfied for t = 0. The rest arguments are also really similar to those 
for the non-relativistic case and we omit them. 

Remark 2.2 For interface problems for the RMHD equations when the interface moves with 
the velocity of plasma particles (e.g., for current-vortex sheets, see Sect. 4) it is useful to have 
the representation 

Aj = VjAo + Gj, (31) 

with 





( 





ej — Vjv'^ 


\ 






-VjV 


O3 






v 







/ 



V4 < 2—— v (E> v 



iv,H) 



v(S)H + H (S)v) 



r r 

+ ^ S^^{v H + H (E) v) -2{v,H){I -v(S)v)^ 



^ ® Cj + Cj (S) H) — — {v (g) Cj + Cj (g) v) 



3 Secondary symmetrization of the RMHD system 

Now our goal is to find a diflFerent symmetrization of the RMHD equations. In some sense it 
will be a secondary symmetrization of system ([8]) which was already symmetric. This secondary 
symmetrization is a relativistic counterpart of the secondary symmetrization of the MHD system 
proposed in [lOJ. In Sect. 4 it will play a crucial role in finding a sufficient stability condition 
for relativistic current-vortex sheets. 

We again consider the linear constant coefficient symmetric system (jl9p - (|2ip / (j22p for the 
perturbation 5V' in the rest frame. For this system we have the standard conserved integral 



d_ 
di 



{A'o{U')6V',6V')dx = 0. 



(32) 



At the same time, we can derive another conserved integral. Indeed, multiplying (j20p and (]2ip 
by 5H' and 6u' respectively, and using ()19p and the divergence constraint div6H' = 0, we get 
the following conserved integral for L2 solutions of the Cauchy problem: 



[ \ir(H',6u')6p-(A'6u',6H')\dx = 0. 
Jr3 [cj J 
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The last integral can be rewritten as 

/ {Co5V',6V')dx = Q, 



(33) 



with 





72" 









"1 







-A' 


) 



Combining (j32p and (j33p . we obtain the new conserved integral 

d 



- / {%{U')5V\5V')dx = ^, 



(34) 



where 



/ 1 ,H'^ 

2 ^—3r 







pa- ci 
X^- A' -XA' 



\ -A^' / / 



and A = X(U') is an arbitrary constant. 

Let us return to the nonlinear setting and consider the nonlinear counterpart 

_ 3 _ 

A'o{U')dtV' + Y,A'jiU')djV' = 



(35) 



of system ([22]) . where V = V\y=Q and U' = U\^=q = {V\y=Q,S) are now the unknowns in the 
rest frame rather than constant vectors. For finding a symmetric system corresponding to the 
conserved integral p4p. where A = X(U') is now a function of U', we multiply (j35p from the left 
by the matrix 



D' = D\U') = %iU'){A',{U')) 



-1 



/ 1 Ah'^ 



XphH' 
\ 





-XA' 



-XA' I I 



One can understand that to get again a symmetric system we should add to the result of 
multiplication the vector R'di\H', with R' = R'{U') = -A(l, 0, 0, 0, i?'). 
Indeed, 



D'A'^ = % = {%) 



f—H' 




\ 




XH'jA' 






JCj 


-XJCJJ 
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That is, 

_ 3 _ 

D'{U')\^Ao{U')dtV' + ^Aj{U')djV'} + R'{U') divif' 



%{U')dtV' + 2i;(C/')5,F' = 0, (36) 



with 



Cj XH'jA' ICj 



and 



Mj = H'J-{H'(^ej + ej(^H'), Cj 



i/j^' 



Note that 



D' > 4^ < 



pa 



2 + ' 



i.e., % > provided that (cf. (fTUl) ') 



p > 0, > 0, < Cs < 1, and A^ < 



To obtain the secondary symmetrization 



MU)dtV + ^%{U)djV = 



(38) 



in the LAB-frame corresponding to ()36p in the rest frame we again use the Lorenz transformation 
(see Sect. 2): 

_ _ 3 _ _ _ 

2lo = j"^r(2l'o + Y,Vk%)j = j'^iCo + Xr{Co + C))J = Ao + XBo = J'^^qJ, (39) 

k=l 



k=l 



Aj + vj XBo + A (vj J^CJ + J"^£j j) , (40) 



where A = A({7), 



k=l 



-V 



pa^ 



{v,H')A' 
V J\f J 



, A/" = (v, F')/ - (//' v + ® i/') 
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f l + X{v,H') T^.H'^ 

7^ V + A TT' 



Co = Co + Ar(£o + C)=T 



-Xv 



T \ 



v + X- 



H' 



[l + X{v,H'))A' IC^-XA' 



\ 



I 



-Xv K-XA' I + XM 

Note that (v, H') = (v, H) and using (|24p all the other values appearing in the formulas above 
can be rewritten in terms of H. System (|38p together with equation (|lip form the symmetric 
system 

3 

" (41) 



%Q{U)dtU + ^'^j{U)djU = 



with 



2ln 



2to 0^ 
, 1, 



21,- 



2t, 

V ^J. 



(42) 



where, if necessary, the symmetric matrices 21^ can be written in a concrete form (in terms of 
dyadic products) like matrices Aa in 

Regarding the hyperbolicity condition 2lo > for system (j^T]) . it is equivalent to the require- 
ment Co > which holds provided that 

p > 0, > 0, < < 1, and < „ ' (43) 

(we drop technical calculations). It is natural that the hyperbolicity condition (|43p for the 
secondary symmetrization (j4ip in the LAB-frame coincides with the hyperbolicity condition 
(1371 ) in the rest frame because 



4 Stability of relativistic current-vortex sheets 

We consider the RMHD equations for t G [0, T] in the unbounded space domain M'^ and assume 
that 

S(t) = {xi=(p(t,x')} (44) 

is a smooth hypersurface in [0,T] x M^, where x' = (a;^,a;^) are tangential coordinates. We 
assume that S(t) is a surface of strong discontinuity for the conservation laws ([2])-([5]), i.e., we 
are interested in solutions of ([S])"® that are smooth on either side of S(t). To be weak solutions 
of ©-([SI) such piecewise smooth solutions should satisfy Rankine-Hugoniot jump conditions. 
Since we are interested only in current- vortex sheets, we do not write down here the RMHD 
Rankine-Hugoniot conditions in a general case which covers, in particular, shock waves. 

As in non-relativistic settings, for relativistic current- vortex sheets we require that the hy- 
persurface S(t) moves with the velocity of plasma particles at the boundary and the magnetic 
field on S(t) is parallel to S(t). With these requirements the Rankine-Hugoniot conditions give 
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the boundary conditions 



dtV = 4, (45) 
[q] = 0, (46) 
i/± = onS(t), (47) 

where [g] = s^Ie — 5 Is denotes the jump of 5, with := g in Q^{t) = {x^ ^ ip{t,x')}; 

= {v^,N), = {H^,N), and N = (1, — (?2y?, — Ssv?)- Again, as non-relativistic settings 
[llj . one can show that (j47p are not real boundary conditions and should be considered as 
restrictions on the initial data 

U^{0,x) = U^{x), xen^iO), (p{0,x') = ipo{x'), x'eM^ (48) 

Remark 4.1 Our assumption that the hypersurface (j44p has the form, of a graph is not a strong 
restriction for special relativity. Regarding the case of general relativity, S should be a compact 
codimension- 1 surface. However, as for shock waves we can follow Majda's arguments 18] for 
extending the results to a compact hypersurface. On the other hand, even for general relativity, 
this assumption is not still a strong restriction because, as in 15], we can consider the interface 
not only locally in time (in the sense of local-in-time existence of piecewise smooth solutions), 
but also locally in space (for compactly supported initial data). 

As in the non-relativistic case jlO^ lllj . the linear stability condition for planar current- 
vortex sheets which should be satisfied at each point of the initial hypersurface S(0) is the basic 
assumption providing the local-in-time existence of smooth (nonplanar) current-vortex sheet 
solutions, i.e., the existence of a solution [U^ ,ip) of the free boundary value problem ([2])- 
dS]), (Il5])-(I18]), where := U in il^{t), and is smooth in n^{t). As for non-relativistic 
current- vortex sheets [10], we formulate the constant coefficient linearized stability problem: 

3 

A^dtU"^ + AfdjU"^ = if X G M^, (49) 

dtip = vf- V2d2ip - vfd-sip, [q] = 0, if x^ = 0, (50) 

U'^{0,x) = U^ix), X £ M|, (/7(0,x') = ifoix'), x' G R^, (51) 

where {U~^ , U~ , (p) is now a perturbation of the constant current- vortex sheet solution (U^ , U" , 0) 
corresponding to a planar relativistic current-vortex sheet with the equation x^ = (all the 
hat values are associated with this solution and without loss of generality we take (p = 0), 
R% = {±x^ > 0, x' £ M^}, Aa = AaiU), [q] = q^^^^^ - q^^i^Q- Since {U+, t/", 0) is assumed to 
be a current-vortex sheet solution, in view of (j45p and (j47p . 

C/± = (p±, 0, vt,vf,0, H^Hf, S). 
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Again, as in non-relativistic settings [lOllllj. we can show that the solution of problem (j49p -([5T 
satisfy the linearized conditions (|T7|) . 



= H2d2f + Hfd^<f if = 0, (52) 



provided that these conditions were satisfied by the initial data (l5T|l . 

The linearized stability of a planar relativistic current- vortex sheet means the well-posedness 
of problem (|49p -(|5i p and is equivalent of the fulfilment of the Kreiss-Lopatinski condition (see, 
e.g., [51 Sj). As was shown in [lOj, for non-relativistic current- vortex sheets the Kreiss-Lopatinski 
condition can be satisfied only in a weak sense, i.e., they can be only neutrally stable. Here we 
do not formally show that the same is true in the relativistic case, but it is natural to expect that 
this is really so, i.e., the uniform Kreiss-Lopatinski condition is never satisfied for relativistic 
current-vortex sheets. Due to technical reasons [10] it is almost impossible to study the linear 
stability of current-vortex sheets by the standard method of normal modes. It is clear that for 
the relativistic case the situation is even worse. Therefore, as in [10], we will use the energy 
method for finding a sufficient neutral stability condition for relativistic current-vortex sheets. 

In the application of the energy method the secondary symmetrization (j4ip found in Sect. 
3 plays the crucial role. The linearization of ()4ip against the constant solution reads: 



3 

^"^dtU^ + ^fdjU^ = if X G M^, (53) 

where 21^ = 2la(C/^). For systems (j53P standard arguments of the energy method give the 
energy integral identity 

^/(t) - ^ [{%iU, U)] l,^^ dx' = 0. (54) 

where 

[{%,U, U)] l,^^= {%tu+, C/+)Ui=o - i^iU-, t/-)Ui=o, = 

I{t) = [ /3+(2i(tC/+,[/+)d3;+ [ l3-{%U-,U-)dx, 

JrI Jri 

and the constants /3+ > and > will be chosen later on. As follows from (|54p . the 
dissipativity condition [(5SiC/, U)] |^i^q< satisfied for an unperturbed fiow [U^, U^) guarantees 
the stability of the corresponding planar current- vortex sheet. Following lllj. we can choose 
such = \{U^) in the secondary symmetrization (|4ip that 

[(^if/,f/)]Li=o=0' (55) 

and the last necessary inequality in the hyperbolicity conditions (j43p for chosen gives us a 
sufficient stability condition. 

Before choosing A^ we have to compute the boundary matrices 21^. Actually, we need not 
these matrices themselves but the quadratic forms (2tf C/^, U^)\^i^q. Let us compute them even 
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for a general basic state {U^{t,x),U {t,x),if{t,x')) satisfying conditions (jl5|) and (|17|) . Then 
the boundary matrices read 



We first calculate the quadratic forms (2l^?7^, C/^)|^i=o ^or = 0, i.e., the quadratic forms 



(^±C/±,C/±)|,i=o- 

Taking into account ([15]) . ([TFl) and representation ([3T]) . we have 

After some algebra we find 



x^=0 



where = (u=^,iV), iV = {1, -820, -dsip), 

and 5(5^)^ is the perturbation of for itx^ > 0. 
Now we consider the case ^ 0: 

u) = {A%U, u) + A± (mt^c/, [/) U^,, 



where 



xi=0 



Omitting (not extremely long) calculations, we obtain 



2(p± + {{H')^, {H't)) ({{H')^, N) + ^t^^ {v^, {H')^) 



That is. 



and for the planar discontinuity with = 

{^fU,U)l,^^ = 2f±5± ( (1 - X^{v^,H^)) vf 



r±({,±)2 

A± 



^2 



a;i=0 



(r±)2 



A^ 



(56) 



xi=0 
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We now choose in (|54p : 



r±(i-A±(?)±,F±)) 

We assume that the hyperbohcity condition (f43l) holds for U = U^, i.e., 21^ > 0. The last 
inequality in (gS]) guarantees that 1 - X'^{v^,H^) > 0, i.e., > 0. Then I{t) > and it 
follows from ()56p that 

[(^i[/,f/)]|^,^o= 2g+[r;i-Ai/i]Ui=o. 

where 

A± = ; (57) 

(r±)2(i-A±(i)±,ii'±)) 



The boundary conditions (jSOll and (l52]l imply 

where -O"^ = {v2,vf), -ff"^ = {H^,H^) and Vtan = (<92,<93). As for non-relativistic current- 
vortex sheets [11] , we choose A^ such that 

- A^-ll] = 0, (58) 

i.e., ()55p holds. This gives us an a priori L2 estimate if the chosen A^ satisfy the last inequality 
in (|43p . At the same time, this inequality for the chosen A^ gives a sufficient stability condition 
for planar relativistic current-vortex sheets. 

Assume that the tangential magnetic fields H^^^ are nonzero and nonparallel to each other: 

H+H^ -H+H^ ^Q. (59) 

For non-relativistic current-vortex sheets the violation of (j59p corresponds to the transition to 
instability. It is natural to expect that the same is true in the relativistic case. It follows from 
(I58P and ([59]) that 

_ |sinv9=F| 



|i?ll±| I sin((^+ -(/?-)! 
where 

mil \H»^ 



In view of ([57|) . 



l + A±(t)±,i/±) 
The last condition in (I43p for A^ gives 



A± = . ^ , with A± = (f±)2A±. 



I A I < m , with m = 



/5±(a±)2 + (B±)2 
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We can show that in terms of the last inequahties read 

1 



A±l< 



I.e., 

l/(f±)^ 



A±I< 



Then, we can finally write down the sufficient stability condition in the same form as for the 
non-relativistic current-vortex sheets (we drop the hats, i.e., write down the condition to be 
satisfied for the initial data in a counterpart of the nonlinear existence theorem from [11]): 



G{U+,U-)>0, (60) 

where 



G{U+, U') = \ sin(v9+ - ip-)\ mm 
but now (in the relativistic case) 



7" 



sm ip~ I I sm if' 



± 



\H\\^\ct[l-{v^f) 



Vp±(a±)2 + (i3±)2+c±|(^±,//±)|' 
Recall that we also assume (cf. ([59])) 

H^H^ - H^H^ > e > 0, (61) 

where e is a fixed constant. 

If in 7^ we formally set = and the specific enthalpies = 1 (this corresponds to the 
non-relativistic limit), then (|60p coincides with the sufficient stability condition for the classical 
current- vortex sheets [lOl [11] . 

Since current-vortex sheets can be only neutrally stable, there appears a loss of derivatives 
phenomenon in a priori estimates. Therefore, the nonlinear existence theorem was proved in 
[llj by Nash-Moser iterations. In the relativistic case, all the arguments towards the proof of 
the local-in-time existence of smooth current-vortex sheet solutions, provided that the stability 
condition ([Un|) , ([HT]) is satisfied at each point of the initial discontinuity, are absolutely the same 
as in the non-relativistic case in |11] . 
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